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ABSTRACT 


Knowledge  of  the  large-scale  vertical  velocity  is 
required  to  predict  the  evolution  of  the  atmospheric 
planetary  boundary  layer  (APBL)  .  Since  naval  operations  are 
often  conducted  in  data  sparse  regions,  single-station 
assessments  of  the  vertical  velocity  are  necessary.  Data  to 
evaluate  such  assessments  were  obtained  from  rawinsondes 
taken  at  San  Nicolas  Island,  California.  Vertical  velocity 
estimates  obtained  by  vertical  integration  of  the  moisture 
budget  equation  (Q-method)  and  by  the  adiabatic  method  were 
used  to  initialize  an  APBL  24-hour  prediction  model.  RMS 
error  statistics  on  predicted  inversion  height,  potential 
temperature,  and  specific  humidity  were  computed  for  the 
forecasts  and  compared  to  the  RMS  errors  of  a  persistence 
forecast-  Calculation  of     the      vertical      velocity   by     the 

Q-method        showed        the        most        promise.  However,  no 

single-station   assessment    method    improved   on   the   persistence 
forecasts. 
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I-       INTRODUCTION 

One  aspect  of  naval  operations  is  battle  group  opera- 
tions in  a  fixed  geographical  area.  This  is  currently 
referred  to  as  MODLOC.  Several  scenarios  lend  themselves  to 
MODLCC  operations:  the  Viet-Nam  era  carrier  based  strike 
operations,  pre-positioning  for  a  projection  of  power,  or 
simply  a  show  of  force.  Training  operations,  as  well,  are 
often  conducted  in  fixed  areas  to  avoid  hazards  to  commer- 
cial shipping  and  air  traffic  with  weapons  firings,  or  for 
privacy   in   operations. 

A  second  aspect  is  the  need  for  self  sufficiency  of  a 
deployed   force.  Any   autonomous     ability   to      make  environ- 

mental prediction  within  the  battle  group  reduces  the 
reguirements  for  shorebased  support,  decreases  delivery 
time  for  perishable  information  and  reduces  the  demand  on 
communication  assets. 

A  third  aspect  is  the  effect  of  the  atmosphere  on  the 
perfcrmance  of  weapons  and  sensors  that  depend  on  electro- 
magnetic and  electro-optical  wave  propagation.  Meteorologi- 
cal factors  which  affect  such  propagation  can  change 
significantly  in  four  to  six  hours.  The  importance  of  APBL 
prediction  to  naval  operations  lies  in: 


•  Prediction   of      the   existence  and  evolution     of  electro- 
magnetic ducts. 

•  Prediction    of  the   inversion    height,        where   it  is    known 
that   optical   propagation   is    degraded  "due   to  turbulence. 

This  thesis  is  an  examination  of  the  accuracy  of  single- 
station  assessments  of  the  large-scale  vertical  velocity, 
and  the  effectiveness  of  these  assessments  for  an  existing 
APBL  prediction  model.  Three  methods  of  computing  the  ver- 
tical velocity  will  be  discussed.  Vertical  velocities  cal- 
culated by  single-station  methods  will  be  used  in  an  APBL 
model  to   ascertain   the  relative    merits  of  these   methods. 
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II.       BACKGROUND 

The  single-station   forecasting   problem    probably  has   been 

of    interest     as   a  military     reguirement  since   the     origin   of 

armed  conflict.      Oliver  and  Oliver   (1945)    put  single-station 

assessments  into  a  naval   perspective: 

During  the  last  few  years,  wartime  conditions  have  made  it 
necessary  for  isolated  combat  units  to  issue  forecasts  in 
regions  where  no  network  of  meteorological  stations  could 
be  available.  Frequently  the  data  from  several  stations 
or  from  reconnaissance  planes  are  available,  but  in  some 
regions  the  forecaster  must  rely  only  on  surface  and 
upper-air  observations  made  at  his  own  station.  This  is 
particularly  true  in  the  case  of  ships  at  sea.  Hence,  it 
is  important  to  develop  proficiency  at  extracting  informa- 
tion with  limited  aerological  data....  He  must  further  be 
able  to  distinguish  from  solutions  which  are  internally 
consistent,  and  those  which  are  inconsistent  and  therefore 
to   be   discarded. 

In  APBL  assessments  the  assumption  of  horizontal  homo- 
geneity (one- dimensionality)  is  often  used.  All  APBL  pre- 
diction models  require,  as  an  externally  known  parameter, 
the  large-scale  vertical  velocity  which  direcly  affects 
changes  in  the  height  of  the  APBL.  Furthermore,  the  height 
of  the  APBL  affects  other  properties  of  the  boundary  layer 
such   as   the   temperature   and  humidity. 

Data  for  this  examination  were  obtained  from  rawinsondes 
taken  at  San  Nicolas  Island  (SNI),  California.  SNI  data  are 
sufficiently  representative  of  a    marine   environment   for  this 
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purpose.  Battalino  et  al.  (19  79)  have  stated  that  condi- 
tions at  Sfll  may  be  marine,  continental,  or  mixed  depending 
on  the  circulation  patterns,  A  decision  to  use  the  SNI  data 
was  made,  in  part,  because  of  the  availability  of  similar 
data  from  Vandenberg  (VND)  Air  Force  Base,  California  and 
Point  Mugu  (PMTC) ,  California,  so  that  examination  of  the 
horizontal  variability  of  the  area  and  divergence  derived 
vertical  velocity  values  could  be  accomplished.  Pig.  1  is 
provided    for   geographic   orientation. 


y«  » 
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Figure   1.      Southern   California  Ccastal   Waters.       (After 
Battalino    et   al-    1979.) 
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A.       THE    PREDICTION    MODEL 

An  APBL  prediction  model  in  which  the  mean  vertical 
motion  is  an  external  parameter  is  used  to  establish  the 
affect  of  the  derived  values  of  vertical  velocity  on  boun- 
dary layer  evolution-  The  model  used  is  based  on  the  one 
presented  by  Stage  and  Businger  (1981)  and  was  coded  for 
operational  use  by  Brower  (1982)  .-  The  predictive  equations 
are: 

dq/dt      =   h-MM1?^)     ♦   We£q)  (1a) 


d9e/dt  =  h-Mw'9e'(0)    +  Rc    "    Rb    +   »e^®e)  (1b) 

dh/dt  =      WD    ♦   We  (1c) 

Where  q   is   the  total   water  content    (liquid   and   vapor)  ,      h   is 
the    inversion   height,    9e    is  the    equivalent    potential  temper- 


ature, w,q,(0)  is  the  surface  flux  of  moisture,  w'9^  (0)  is 
the  surface  buoyancy  flux,  A  is  the  difference  between  the 
variable  value  immediately  above  the  APBL  and  its  well  mixed 
value  in  the  APBL,  Rc  is  the  rate  of  radiative  heat  gain  per 
unit  area  at  the  cloud  base,  and  Rb  is  the  rate  of  radiative 
heat  loss  per  unit  area  near  the  cloud  top.  We  is  the 
entrainment  velocity  and  WD  is  the  large-scale  vertical 
velocity    at  the   inversion   height. 
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The  model  atmosphere  is  two  layered,  composed  of  a  tur- 
bulent lower  layer  in  which  equivalent  potential  tempera- 
ture, specific  humidity  and  wind  speed  are  constant  with 
height,  and  a  stable  upper  atmosphere  with  constant  specific 
humidity    and   buoyancy  gradients.  The  layers   are  separated 

at  the  inversion  height  by  a  zero-order  discontinuity 
(jump).  Fig.  2  depicts  how  a  sample  sounding  is  character- 
ized  in  the   model  atmosphere. 

Surface  fluxes  of  moisture  and  bouyancy  are  diagnosed 
from  bulk   formulae    (Lilly    1968).  The  liquid   water   profile 

is  computed  adiabatically  and  reduced  by  about  30  percent  to 
agree  with  empirical  results  (Fairall  et  al.  1981.)  Long 
wave  radiative  cooling  is  treated  using  the  Stephan-Boltzman 
law,  and  short  wave  radiative  warming  is  treated  using  the 
delta-Eddington    approximation    (Fairall  et   al.      1981). 

The  large-scale  vertical  velocity,  WD  ,  is  treated  as  an 
external  parameter,  and  is  constant  over  the  model  run 
(2«*-hr)-  The  system  is  closed  with  the  specification  of  the 
entrainment  velocity  which  is  accomplished  by  the  assumption 
that  dissipation  of  turbulent  kinetic  energy  (TKE)  is  a 
fixed  fraction  of  the  production  of  TKE  (Stage  and  Businger 
1981).      A    30-minute  time  step  is    used. 
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Figure  2.   Sample  sounding  in  the  model  atmosphere.  Solid 

line  represents  data,  dashed  line  represents  the 
model  atmosphere. 
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III.       MET  HODS 

At  least  three  methods  which  can  be  used  to  compute  the 
large-scale  vertical  velocity  are: 

•  The   kinematic  method. 

•  The   adiafcatic  method. 

•  Integration  of  the    moisture    budget   equation. 

The  kinematic  method  requires  wind  sounding  data  from  at 
least  three     appropriately  spaced  stations.  The  adiabatic 

methcd  and  the  integration  of  the  moisture  budget  equation 
are   truly   single-station  assessments. 

A.       THE    KINEMATIC    METHOD 

The  kinematic  method  can  be  used  to  compute  the  vertical 
velocity  when  nearly  simultaneous  (within  two  hours)  sound- 
ings are  available  at  three  or  more  appropriately  spaced 
stations  such  as  VND,  PMTC,  and  SMI  (Fig.  1).  Following 
Saucier    (1955),    the   vertical  velocity   can    be   derived: 

dw/c-z     =   -    (du/dx      ♦        dv/dy)  ,  (2) 

where  u,  v  and  w  are  the  zonal,  meridional  and  vertical  wind 
components.  The  horizontal  wind  components  can  be  linearly 
approximated   by   a   Taylor  series    expansion   as: 

u   =  u.    ♦    (du/dx)    x    +    (du/dy)     ?*  (3a) 

v   =  v0    +    (dv/dx)    x   +    (dv/dy)    y.  (3b) 
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Eg.  (3)  suggests  that  the  horizontal  derivatives  may  be 
obtained  from  a  regression  of  the  horizontal  wind  components 
onto  the  horizontal  plane  (Fairall  et  al.  1981)  .  Decomposi- 
tion of  the  wind  vector  from  the  sounding  data  at  the  three 
stations  (SNI,  VND,  PMTC)  into  zonal  and  meridional  compo- 
nents is  required.  Vertical  profiles  of  the  horizontal  wind 
components  were  constructed  at  a  200  meter  vertical  incre- 
ment by  linear  interpolation  from  the  irregularly  spaced 
data  levels  to  those  levels  that  are  an  integer  multiple  of 
200  meters,  Multiple  linear  regression  of  each  component  of 
the  horizontal  wind  onto  the  horizontal  plane  is  performed 
as  fellows: 

u  =  a  ♦  bx  +  cyr  (4a) 

v  -  d  ♦  ex  ♦  fy,  (4b) 

where  a  through  f  are  the  regression  coefficients. 

Comparison  of  Eg.  (3)  with  Eq.  (4)  indicates  that  the 
coefficients  b  and  f  give  the  zonal  derivative  of  zonal  wind 
and  the  meridional  derivative  of  meridional  wind  at  the  cen- 
ter cf  a  triangle  formed  by  the  three  locations  and  at  each 
vertical  level  where  the  regression  is  done. 

Vertical  integration  of  Eq.  (2)  is  performed  to  give  a 
vertical  profile  of  vertical  velocity  to  three  kilometres: 

3  km 
w(z)  -  -   £((b)   ♦  (f)  )  Az  f  (5) 
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where  w(z=0)  *  o  and  Az  =  200m. 

B.   THE  ADIAEATIC  METHOD 

For  isentropic  motions  there  is  no  time  change  of  poten- 
tial temperature  following  the  motion.  Based  on  this 
assumption  Saucier  (1955)  specifies  the  adiabatic  method  for 
estimating  vertical  velocity  as: 

(g9/dt)  +  V  (V9) 

v  ~      -    ,      (6) 

oe/az) 

where  8  is  the  potential  temperature.  Saucier  further 
states  that  the  local  temperature  change  due  to  advection 
may  be  on  the  order  of  that  due  to  vertical  motion.  In  this 
study,  calculations  of  vertical  velocity  are  made  both 
including  and  neglecting  the  advective  term.  Description  of 
the  advection  calculation  is  left  to  a  subsequent  section. 

The  vertical  gradient  is  determined  by  a  linear  fit  to 
the  lapse  rate  of  potential  temperature  above  the  inversion. 
The  local  time  change  of  temperature  is  obtained  at  fixed 
heights  (z  =  1.0,  1.25,  1.5  km)  from  sequential  soundings  at 
the  same  location.  These  heights  are  above  the  inversion. 
It  is  assumed  that  the  vertical  velocity  at  the  surface  is 
zero.  The  vertical  velocity  at  the  inversion  is  taken  as 
the  average  of  the  vertical  velocity  at  the   fixed  heights 
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linearly   scaled   to  the   inversion      height.        Thus,      a  single- 
station   assessment  of  the   vertical  velocity   is  obtained. 

C.       INTEGRATION    OF   THE    MOISTURE    BUDGET 


Integration  of  the  moisture  budget  equation  to  compute 
the  vertical  velocity  was  performed  by  Lenschow  (1973). 
This  method  will  hereafter  be  referred  to  as  the  Q-method, 
and  is  based  on  the  assumption  of  well-mixed  specific  humid- 
ity in  the  boundary  layer,  and  that  changes  occur  due  to 
fluxes  at   the  sea  surface   and  inversion  only. 

Vertical  integration  of  the  moisture  budget  equation 
through  the  boundary  layer  and  across  the  interface  between 
boundary  layer  and  the  stable  atmosphere  above  yields  the 
following  expression  for  vertical  velocity,  w,  at  the  level 
of    the  inversion,   h: 


A(q+1)  (dh/dt)     -    h(d(q+l)/dt)     ♦    w«q«(0) 

w   =  ,         (7) 

Mq+D    ♦    (0h/2) 


where  q  and  1  are  vapor  and  liquid  water  contents,  w*qf  is 
the  moisture  flux,  /?  is  the  vertical  gradient  of  (q  +  1)  in 
the  mixed  layer,  and  A(q*l)  is  the  difference  between  total 
moisture  in  the  mixed  layer  and  total  moisture  immediately 
above  the  mixed  layer.  With  no  liquid  water,  z  «h,  and 
neglecting  moisture  advection  Eg.  (7)  is  simplified  to: 
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Aq(dh/dt)    -    h(dq/dt)    ♦    w«g»(0) 

w   =  .  (8) 

Aq 


Lilly  (1968)  assumes  that  (5  is  zero.  Lenschow  (1973)  found 
£h/2  to  be  small  and  certainly  much  less  than  the  moisture 
jump. 

Eq.  (8)  provides  another  independent  assessment  of  the 
vertical  velocity  at  the  inversion.  The  time  derivatives 
(dh/dt  and  dq/dt)  are  determined  from  sequential  soundings. 
The  moisture  jump  and  inversion  height  are  taken  as  averages 
for  the  two  soundings,  and  the  surface  moisture  flux  is  det- 
ermined  by   the   bulk   method,   i.e.. 


Wq'(O)     =    ChV1Q(qs     -    q1Q)  (9) 

The  exchange  coefficient,  Ch,  is  computed  from  the  wind 
speed  dependent  relation  of  Large  and  Pond  (1980) .  The  sur- 
face specific  humidity,  qs,  is  determined  via  Tetens*  for- 
mula  for   the    vapor   pressure    (Buck   1981). 

D.        DETERMINATION    OF    ADVECTION    FROM    THE   THERMAL    WIND 

Under  the  geostrophic  assumption,  the  horizontal  temper- 
ature gradient  is  directly  proportional  to  the  vertical 
shear  of   geostrophic   wind-      Following   Byers    (1974): 

v  =      f"Mc-*/dx)p  ;      u   =  -f-Mc-*/oT)p  .  (10) 
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where  v  and  u  are  the  meridional  and  zonal  components  of  the 
geostrophic  wind,  f  is  the  Coriolis  parameter,  <t>  is  the 
geopctential  and  the  horizontal  derivativss  are  on  the  con- 
stant pressure  surface.  Vertical  differentiation  of 
Eq.  (10)  and  substitution  from  the  hydrostatic  equation 
yields: 

dv/dp  =   f-»  (-a  )    ,  (11a) 

du/3p  -*-*(«),  <11b> 

where      a    is  the     specific   volume.  Substitution   from      the 

equation  of  state  and  rearranging  terms,  the  components  of 
the    temperature   gradient  are: 

(c-T/dDp      =      -(fp/R   )  (c-v/dp)     ,  (12a) 

(dT/dy)p      =         (fp/R   ) (du/dp)    ,  (12b) 

where  R   is   the  gas  constant   for    dry  air. 

With  a  wind  sounding,  assumed  to  be  geostrophic  above 
the  boundary  layer,  the  horizontal  temperature  gradient  can 
be  specified  and  a  single-station  assessment  of  the  thermal 
advection  term  in  Eq.  (6)  is  obtained:  -  (V- VT)  =- (V- V9)  ,  at 
a   constant    level. 

This  calculation  is  performed  by  subtracting  the  wind  at 
the  first  level  above  the  inversion  from  the  wind  at  an 
upper   level.  The   advecting   wind   is     the    mean   wind      in   the 

layer. 
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E.       SINGLE-STATION    ASSESSMENT   SCORES 

The  separate  vertical  velocities  computed  from  single- 
station  assessments  by  the  adiabatic  method  and  the  Q-method 
were  both  used  as  the  large-scale  vertical  velocity  in  the 
APBL    model. 

Eoth  methods  required  use  of  finite  differences  in  time 
to  approximate  time  derivatives.  The  time  increment  used 
was  the  time  difference  between  subsequent  soundings  at  SNI; 
it  varied  from  four  to  eighteen  hours.  To  establish  the 
vertical  velocity  for  a  model  run,  the  time  weighted  average 
vertical  velocity  was  used,  as  computed  from  the  SNI  data 
for    the   24    hour    period   following    model    initialization. 

The  model  requires  initial  values  of  inversion  height, 
and  initial  profiles  of  potential  temperature  and  specific 
humidity.  Sea  surface  temperature,  required  for  computing 
surface  moisture  flux,  was  constant  through  the  model  day. 
Wind  speed,  required  for  computing  surface  moisture  and 
bouyancy  fluxes,  was  input  as  a  linearly  changing  function 
of    time. 

Three  model  runs  were  made  on  each  set  of  initial  condi- 
tions: one  with  vertical  velocity  by  the  Q-method,  a  second 
with   vertical  velocity   by   the      adiabatic   method   and  the   last 
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with  vertical  velocity  set  to  zero.  The  advective  term  of 
Eg-  (6)  it  was  found  to  be  much  less  than  the  local  time 
change  of   potential  temperature.      Thus,    it    was   neglected. 

Foot  Mean  Square  (RMS)  error  was  computed  for  predicted 
inversion  height,  hr  potential  temperature,  0,  and  specific 
humidity,    q. 

_J1/2 

RMS  error  = 


i|§.,   -*ff 


where  N  is  the  number  of  observations,  xt  is  the  observed 
value  cf  a  variable,  and  y,  is  the  predicted  value  of  a 
variable.  Observed  values  were  obtained  from  rawinsonde 
data  available  within  the  24  hours  following  the  time  of  the 
initial  data.  Typically,  two  or  three  soundings  were  avail- 
able for  each  24-hour  period  following  initialization. 
Predicted  values  were  taken  from  the  model  output  for  the 
time  co  rr.es  pen  ding  to  that  of  the  observed  values  (±  15 
minutes,    due   to   the   30   minute  model   time  step). 
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IV.   SYNOPTIC  SITUATION 

Two  periods  (cases)  were  considered.  Case  I  was  from  31 
August  to  4  September  1981  and  consisted  of  the  analyses  of 
nine  model  predictions  verified  against  twenty  observations. 
Case  II  was  from  15  to  18  September  1981  and  consisted  of 
the  analyses  of  eight  predictions  verified  against  seventeen 
observations. 

A.   CASE  I  (31  AUG  -  4  SEPT  1981) 

Throughout  the  period  the  eastern  North  Pacific  Ocean 
was  under  a  high  pressure  system  that  was  moving  slowly 
eastward.  Southern  California  was  dominated  by  a  thermal 
trough,  leading  to  weak  offshore  flow  at  the  surface,  while 
California  was  to  the  south  of  the  strongest  westerlies  at 
upper  levels.  Pigs.  3  through  6  show  the  National  Meteoro- 
logical Center  (NMC)  surface  and  500  millibar  (mb)  analyses 
for  Case  I. 

Selected  prints  of  the  GOES  WEST  Infra-red  (IR)  imagery 
indicate  low  broken  to  overcast  cloudiness  in  the  vicinity 
of  SNI  (Fig.  7).  Surface  observations  at  SNI  indicated 
overcast  conditions  with  fog  at   night,   and   scattered  low 
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cloudiness  with  haze  persisting  after  the  late  morning 
burn-off   of  the   fog. 

A  Pacific  high  pressure  system  dominated  the  eastern 
ocean  area  west  of  California  on  31  August  (Fig.  3).  A  low 
pressure  system  was  located  in  the  Gulf  of  Alaska  with  its 
associated  weak  frontal  system.  Tropical  Storm  Irwin  was 
dissipating  southwest  of  the  Baja  Peninsula.  A  thermal 
trough  dominated  the  central  California  valley.  Very  weak 
westerly  flow  aloft  was  predominant  over  southern  California 
and    most   of  the   eastern   Pacific   Ocean. 

After  31  August  the  Pacific  high  intensified,  while  the 
low  pressure  system  in  the  Gulf  of  Alaska  moved  into  the 
northwest  United  States  and  the  assosiated  cold  front  dissi- 
pated. The  thermal  trough  over  inland  California  moved 
slightly  eastward  and  dominated  the  high  desert  region  of 
Arizona  and  southeastern  California  (Figs.  4  through  6). 
Flow  aloft  remained  very  weak  over  southern  California  and 
the  low  level  flow  in  the  vicinity  of  SNI  was  very  weak  (3 
to   5    kt)    and   predominantly   offshore. 


25 


00 


60 

< 


CD 

o 

o 


CO 
0) 
CO 

c 
.a 

e 

o 
o 

Lf) 

c 

CD 

o 
03 

& 

en 


CD 

60 

■H 


27 


GO 


P. 
0) 
00 

CD 

S 
CD 

O 

o 

CN 


<u 

0> 

CO 

to 
td 

£ 

CO 


0) 

W) 
•H 
Cm 


29 


29 


1245  GMT  31  AUG  81 


1315  GMT  1  SEP  81 


1245  GMT  2  SEP  81 


1315  GMT  3  SEP  81 


Figure  7.       GOES    HEST    13,    CASE    I. 
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B.   CASE  II  (15-18  SEPT  1981) 

Throughout  the  period  the  polar  front  was  located  in  the 
eastern  part  of  the  North  Pacific  Ocean  and  Gulf  of  Alaska. 
A  high  pressure  system  dominated  the  area  southeast  of  the 
polar  front.  A  thermal  trough  was  located  over  the  inland 
California  valley  and  through  much  of  the  southwestern 
United  States  and  Mexico.  Figs.  9  through  13  show  the  NMC 
surface  and  500  mb  analyses  for  Case  II. 

Satellite  imagery  indicates  considerable  clear  area  in 
the  vicinity  of  SNI.  Surface  observations  at  SNI  reported 
scattered  low  and  high  clouds  with  occasional  reduced  visi- 
bility due  to  fog  and  haze.  Selected  prints  of  the  GOES 
WEST  IR  imagery  are  shown  in  Fig .  14  . 

Cn  1 4  September  high  pressure  prevailed  over  the  eastern 
Pacific  Ocean.  The  thermal  trough  dominated  the  southwest- 
ern United  States  and  the  Mexican  Plateau.  A  low  pressure 
system  with  a  weak  frontal  system  existed  in  the  Gulf  of 
Alaska.  Weak  unorganized  flow  existed  south  of  the  strong- 
est westerlies  in  a  large  region  of  very  flat  height  gradi- 
ents over  most  of  California  and  the  eastern  Pacific  Ocean 
(Fig-  9)  . 
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Ey  1 5  September  a  large  amplitude  ridge  aloft  had  built 
over  the  western  United  States  and  Canada  (Fig.  10)  .  The 
thermal  trough  weakened  in  central  California  as  the  ridge 
axis  in  the  eastern  Pacific  Ocean  migrated  eastward.  Low 
level  flow  in  southern  California  was  offshore  and  stronger 
(12    to  22   kt)    than  in  Case  I. 

C.       DIFFERENCES    BETWEEN    CASES    I    AND   II 

Comparison  of   the   mixed  layer  and  overlying   data   clearly 

indicates   that     different   synoptic     situations   existed        for 

the    two   cases.        These   differences  are   summarized   in   Table   I 

and    illustrated    in   Fig.    8    : 

TABLE    I 
Comparison   of   the  two  cases 

Case  I  Case   II 

mean    h  604    m                                262    m 

range    of  h  511    to   66  7   m  145    to    357   m 

mean   0  16.2   C                             18. 1    C 

mean  A9  15.0   C                             14. 1    C 

mean   q  8.3    g/k g                         9.8    g/kg 

mean  Aq  -6.0   g/kg  -4.4    g/kg 

mean   d9/dz  above   h  4.2   C/k m                        2.8   C/kra 

mean  dg/dz  above   h  -.35   g/(kg-km)  -.83   g/(kg-km) 

On  the   average,        inversion    height   was   nearly     2.5    times 

greater   in   Case    I   than   in   Case  II.      Case  II    had   a   warmer  and 

more   moist   boundary     layer    with    smaller  magnitudes      of    jumps 
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in  bcth  potential  temperature  and  specific  humidity.  The 
vertical  gradient  of  potential  temperature  (stability)  above 
the  inversion  is  1.5  times  greater  in  Case  I  than  Case  II, 
and  the  vertical  moisture  gradient  is  nearly  2.5  times 
larger  in  Case  I  than  Case  II. 
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Figure  8.      Average   Model  Atmosphere:      Case    I    vs.    Case    II 
Solid:    Case  I,    dash:      Case  II 
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Figure  14.      GOES   tfEST  IS,    Case  II. 
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V.       RESgiTS 

Direct  measurement  procedures  do  not  exist  to  obtain 
vertical  velocity  values.  Therefore  an  indirect  evaluation 
procedure  was  nesessary.  This  was  accomplished  by  using  the 
RMS  error  (section  III,E.)  for  model  predictions  of  inver- 
sion height,  potential  temperature  and  specific  humidity, 
with  the  large-scale  vertical  velocity  specified  by  the 
adiatatic  method,  Q- method,  vertical  velocity  equal  to  zero, 
and  a  persistence  forecast.  It  is  noted  that  the  Q-method 
is  based  on  measured  values  of  the  prediction  variables 
(Eg-  (1)  )  .  Therefore  the  model  is  expected  to  perform  well 
with   vertical  velocities   computed   in   this    way. 

In  both  cases  the  inversion  persisted  for  several  days. 
Eq.  (1c)  illustrates  that  H  acts  to  keep  a  lid  on  the  boun- 
dary layer  only  when  subsidence  (W  <0)  is  occurring  since  W 
provides  only  for  increasing  the  inversion  height,  i.e., 
"detrain ment"  is  not  allowed.  Then  from  a  pragmatic  view, 
negative  vertical  velocities  (subsidence)  are  sought  while 
positive      ones  are     not.  Large-scale     subsidence  is      also 

expected,    frcm  a   synoptic   viewpoint   in   both    cases,    since   SNI 
is    under   the   summertime    subtropical   ridge. 
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Table  II  gives  the  averages  of  the  vertical  velocities 
for  each  case  calculated  by  the  Q-method,  the  adiabatic 
method,  and  the  kinematic  method.  The  standard  deviation  is 
given  in  parentheses-  Figs.  15  and  16  show  the  frequency  of 
occurrence  of  calculated  vertical  velocity  for  Case  I  and 
Case  II,  respectively.  Assuming  that  subsidence  existed  at 
the  inversion  height,  the  adiabatic  and  kinematic  methods 
are  unable  to  diagnose  subsidence.  The  Q-method  consis- 
tently gives  subsidence  where  subsidence  is  expected.  Dis- 
cussion of  possible  sources  of  errors  in  these  methods  is 
left  to   a   subsequent  section. 

TABLE    II 
Vertical  velocity    (cm/sec)    comparisons 

Q-method  Adiabatic  method  Kinematic    method 

Case   I         -.26     (.  14)  0.0    (.22)  1.2    (1.7) 

Case   H      -1.2    (1-1)  -1    (.24)  .7    (1-3) 


A.   DISCUSSION  OF  CASE  I 

The  EMS  errors  for  Case  I  appear  in  Table  III  .  The 
average  time  from  initialization  to  verification  for  this 
case  was  15.6  hours.    Persistence  forecasts  yield  the  least 
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RMS  error  for  all  variables  (h,  9,  and  q)  .  When  the  model 
is  initialized  with  the  large-scale  vertical  velocity  com- 
puted by  the  adiabatic  method  or  when  the  vertical  velocity 
is  assumed  to  be  zero,  very  poor  results  are  obtained  in  the 
inversion   height   prediction.  Vertical  velocities  computed 

by  the  Q-method  result  in  better  prediction  of  inversion 
height  than  the   adiabatic   method. 

TABLE  III 
RHS   Error    (Case   I) 

.6  g/kg 

.8  g/kg 

.9  g/jcg 

.7  g/kg 

The  large  errors  in  the  inversion  height  for  the  W=0 
prediction  occur  because  the  boundary  layer  is  constrained 
to  growth  by  entrainment.  The  poor  results  of  the  adiabatic 
method  are  also  to  be  expected  because  of  two  basic  problems 
with  this  method.  First,  examination  of  Eq.  (6)  reveals 
that  the  lapse  rate  is  assumed  to  be  constant  and  the  numer- 
ator is  approximated  by  a  finite  difference  in  time  over  six 
to  twelve  hours.  An  uncertainty  of  ±  one  degree  in  poten- 
tial temperature  can  lead  to  an  uncertainty  in  vertical  vel- 
ocity of  ±.5  (cm/sec)  ,  assuming  nominal  values  of  the  free 
atmosphere   lapse      rate    (2.5  to      5   C/Km)  .         Since      the    actual 
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vertical  velocity  may  be  positive  or  negative  and  often  on 
the  order  .5  (cm/sec),  the  adiabatic  method  is  not  adequate 
to  specify  either  the  magnitude  or  the  sign  of  the  vertical 
velocity-  Second,  one  must  choose  an  appropriate  time 
interval  over  which  the  finite  difference  applies.  In  prac- 
tice the  time  interval  was  that  of  the  rawinsonde  schedule; 
however,  the  soundings  must  be  scheduled  suffiently  close  in 
time  to  reveal  changes  in  sign  of  the  tendencies  of  the 
boundary  layer  variables. 

The  Q-method  also  suffers  from  the  disadvantages  of 
finite  difference  in  the  specific  humidity  (dq/dt)  and  the 
inversion  height  (dh/dt) .  However,  the  moisture  flux  term 
has  no  such  problem.  In  addition,  specification  of  the 
moisture  and  inversion  height  variation  with  time  required 
for  these  calculations  does  not  require  profiling  devices 
such  as  radiosondes,  but  may  be  obtained  from  surface 
instrumentation. 

B.   DISCUSSION  OF  CASE  II 

Average  time  from  initialization  to  verification  was 
14.6  hours.   Table  IV  summarizes  the  results. 

None  of  the  predictions  did  particularly  well  with 
respect  to  RMS  error  of  inversion  height.   The  same  ordering 
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TABLE  IV 
RMS  Error  (Case  II) 

Persistence  1 05  m  1.5  C  . 7g/kg 

H=0  3  75    n  1-5  C  1-4g/kg 

Q-method  142   i  1-0   C  2.4g/kg 

Adiabatic   method  5  27   m  1.8  C  1.3g/kg 

of  RMS  errors  for  inversion  height  prediction  occurs  with 
the  vertical  velocity  specification  methods  in  this  case.  A 
major  difference  in  these  results  is  that  the  model  pre- 
diction of  specific  humidity  with  vertical  velocity  by  the 
Q-method  exhibits  a  21f  error,  much  greater  than  any  other. 
This  is  to  be  expected  since  the  average  inversion  height 
was  262m,  over  40*  lower  than  in  Case  I,  and  in  Eq.  (1a)  the 
specific  humidity  tendency  is  inversely  proportional  to 
inversion  height.  Hence,  low  initial  inversion  heights,  as 
was  the  case,  coupled  with  predictions  of  inversion  heights, 
which  are  too  low,  contribute  to  the  over  prediction  of  spe- 
cific humidity.  The  Q-method  yields  the  least  RMS  error  in 
potential   temperature  predictions. 

C.       CISCOSSICN   OF   THE    KINEMATIC    METHOD 

Vertical  velocities  calculated  by  the  kinematic  method 
were  positive  more  often  than  not  as  indicated  in  Figs.  15 
and  16  .  They  were  much  larger  in  magnitude  than  those  cal- 
culated     by   either      the   adiabatic     method      or   the     Q-method. 
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Therefore  application  of  these  to  the  model  would  have  pro- 
duced larger  RMS  errors  than  any  presented.  Two  reasons  for 
the    poor   results  are: 

•  Improper  specification     of    the     bottom   boundary      condi- 
tion. 

•  Small  sample   size   for  the  regression   approach, 

Cver  the  ocean  the  bottom  boundary  condition  for  Eq.  (2) 
is  that  the  vertical  velocity  equals  zero  at  the  sea  sur- 
face. This  is  a  good  assumption  since  the  surface  is  flat 
with  respect  to  a  local  earth  tangent  plane.  However,  since 
SNI,  PMTC,  and  VND  are  coastal  stations,  the  boundary  condi- 
tion for  Eq.  (2)  requires  specification  of  the  terrain 
(Dutton    1976). 

D.        CONCLUSION 

Since  both  Case  I  and  Case  II  represent  small  samples  of 
observations  the  error  statistics  have  little  stability. 
However,  they  do  suggest  that  vertical  velocity  calculated 
by  the  Q-method  has  the  most  merit  as  a  single-station 
assessment.  Clearly,  the  assumption  of  no  vertical  motion 
or  vertical  motion  derived  from  the  adiabatic  method  has 
questionable  value  in  APBL  assessments.  Fig.  17  shows  a 
model     prediction   with      vertical      velocity      computed   by      the 
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Q-method.  Fig.  18  is  the  same  model  run  but  with  upward 
vertical  velocity  computed  by  the  adiabatic  method.  Compar- 
ison of  these  figures  shows  the  good  agreement  in  g  pre- 
diction indicated  in  Table  III  .  On  the  other  hand,  the 
superior  result  of  the  Q-method  is  well  illustrated  by  com- 
parison of  the  evolution  of  the  inversion  height.  The 
importance  of  diagnosing  subsidence  when  it  exists  is  evi- 
dent  by   this   comparison. 

The  assumption  that  SNI  represents  a  marine  environment 
is  a  possible  source  of  error  in  the  Q-method.  Alsor  sea 
surface  temperature  has  an  indirect  effect  on  the  calcula- 
tion of  vertical  velocity  through  Eg.  (9).  The  sea  surface 
temperatures  used  were  obtained  from  the  Sea  Surface  Thermal 
Analyses  of  the  National  Oceanographic  and  Atmospheric 
Administration,    rather   than  being   in-situ  measurements. 

Additional  research  is  needed  in  specification  of  the 
large-  scale  forcing  of  the  APBL.  Vertical  velocity  by  the 
Q-method  appears  promising,  however  larger  samples  must  be 
obtained  to  gain  statistical  stability.  This  need  not  be  an 
expensive  venture.  One  only  needs  continuous  measurements 
of  sea  surface  temperature,  air  temperature,  humidity  and 
wind.        That      can   be      done   on     a    near-continuous      basis   with 
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currently  available   instrumentation  that   does  not   place 
great  demands  on  personnel. 
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